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Abstract 

We test the fluctuation theorem from measurements in turbulent flows. We study the time 
fluctuations of the force acting on an obstacle, and we consider two experimental situations: the 
case of a von Karman swirling flow between counter-rotating disks (VK) and the case of a wind 
tunnel jet. We first study the symmetries implied by the Gallavotti-Cohen fluctuation theorem 
(FT) on the probability density distributions of the force fluctuations; we then test the Sinai scaling. 
We observe that in both experiments the symmetries implied by the FT are well verified, whereas 
the Sinai scaling is established, as expected, only for long times. 
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I. INTRODUCTION 



The fluctuations of global quantities in out-of-equilibrium systems is a subject of current 
interest which presents many unsolved problems. One of these problems is the prediction of 
the shape of the probabUity deus.ty funetion (PDF) of the variable uuder study; the PDF 

is of course not necessarily Gaussian Among_several approaches used to study these 
problems, the theory of Gallavotti and Cohen 

mm 

, based on a chaotic hypothesis, 
leads to interesting predictions. The main prediction of the Gallavotti- Cohen fluctuation 
theorem (FT) concerns the probability density function of a variable related to the phase 
space contraction rate of the nonequilibrium system under study, this variable being a current 
j such as the flux of heat, or of momentum, or energy, etc. The mean of j on a time interval 
r is defined as 

1 /"^^^ 

Jr = - mdt' . (1) 

T Jt 

One is then interested in the the probability distribution 7!'t{Y) of the variable Y = Jr/Joo, 
where Joo = limT-^oo Jt is the stationary average current studied in the system. If r is larger 
than a characteristic time of the system, then the chaotic hypothesis predicts that tTtCY) 
verifies 

where a is proportional to the phase space contraction rate. It is important to stress that the 
above equation is valid for all values of Y. A related expression for the probability density 
TTriY) has been established by Sinai under quite general assumptions: nr{Y) = e~^^^ . 
In this formulation the function C(^) is a limiting function, for r oo, 

^^^^ ^ ln[7r.(y)]-ln[A.] ^ 

independent of r However, the FT (with the addition of the chaotic hypothesis) is 

giving more information: the function (C(^) — C(~^)) is linear in Y and its derivative is 
related to the contraction rate in phase space. In this sense, the FT takes into account the 
underlying dynamics taking place in the system. 

However, in the original derivation of the fiuctuation theorem, the proof requires several 
very restrictive hypothesis, among which the time reversibility of the system. Therefore it 
is important to check the symmetry predicted by eqI2]in more realistic systems in order to 
probe the universality and applicability in practical cases. Eq|2]has been tested with success 



on a numerical simulation of a rather artificial system ItII, and it is also quite well verified for 

n Q 

the heat flux in a chain of coupled non-linear oscillators j8|, |9| and in the energy dissipation 
rate in the shell model of turbulence 19|1. 

n 

The situation is less clear for local variables. Gallavotti [4] suggested that the predictions 
of the FT could be extended to "the average of local observables" . This has prompted us 
to check this hypothesis using local measurements in a turbulent flow. We started with 
preliminary tests using the data of turbulent Rayleigh-Benard convection These tests 
did show that the fluctuations of a quantity proportional to the heat flux verify both eq|21 
and eqEl In this article we consider the case of the pressure fluctuations in a turbulent 
wind. Specifically, we consider the statistics of the force applied on an object. We find that 
its fluctuations verify the Gallavotti-Cohen predictions. 



II. EXPERIMENTAL APPARATUS 



The experimental set-up which is shown in fig^ consists of two coaxial, counter rotating 
disks. This experimental set-up, known as the von Karman geometry, produces an intense 
turbulence in a compact region of space. It has been previously used to study the fluctuations 
of the power injected to maintain the turbulent flow at a fixed integral Reynolds number — 
more details on the experimental apparatus can be found in ref. 13[. 
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FIG. 1: von Karman experimental set-up. The discs with radius 2R = 20 cm are set H = 20 cm 
apart. They are counter rotated at the fixed frequency $7 = 40 Hz. 
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The distance between the disks is 20 cm. As illustrated in figH] a square plate, of side 
4 cm and 1 mm thick, is inserted between the two disks at a distance of 6 cm from one of 
the two disks, with its center at 6 cm from the disks rotation axis. The plate is mounted 
on an extreme of a 10 cm long lever. The other extreme of the lever is fixed on a strain 
gage (ENTRAN ELFS-T3M-10N force transducer). The lever axis, which is at distance 
of 1 cm from the strain gage, is mounted on very precise ball bearings in order to reduce 
friction. With this system the minimum detectable force exerted on the plate is 1 mN, with 
a response time of 10^^ s. 

The quantity under measurement is the integral of the pressure over the area S of the 
object: 



where p is the fluid density and diUj the velocity gradient tensor. It thus instantaneoulsy 
relates to the strain and vorticity fluctuations in the flow. Note that the presence of the 
object introduces a no-slip condition on its surface which changes the velocity distribution; 
the fluctuations of force on the objet are related — but not strictly equal — to the flux of 
momentum in the undisturbed flow in the absence of the obstacle. The mean value Fq of 
the force on the object is set by the pressure head {p ~ pU"^), where U is the mean velocity 
at the location of the object. The time fluctuations of the force are related to the small 
scale turbulent fluctuations, as indicated by equation (0) and shown in several studies on 
the fluctuations of the local pressure [l^ Q • 

We analyze a set of data recorded at a disk rotation frequency of 40 Hz, corresponding 
to a Reynolds number of about Re = 2 10^. The signal of the instantaneous force F{t) 
is filtered at 70 Hz and sampled at 160 Hz with 21bits resolution Agilent 1430 digitizer 
(tq = 1/160 = 6.25 ms is the sampling time). The signal is recorded for about 2 10^ integral 
times. In these conditions the mean force Fq on the plate is 0.12 N. 

The variable of interest is y = F{t)/FQ. Its typical time evolution is plotted in figI21 in 
(a) is shown the 'standard' fluctuations of the signal in a large time-window while in (b) we 
have selected a time interval during which a negative fluctuation of large amplitude occurs. 
One observes that y = F{t)/FQ fluctuates mainly about the mean value 1, but also that 




(4) 



The pressure itself p(r, t) is the solution of the Poisson equation 



Ap(r,t) = 2pdiUj{r,t)djUi{r,t) 



(5) 
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FIG. 2: (a) Force of the wind on the measuring surface as a function of time, (b) Expanded view 
of the force signal — note the strong negative peaks, (c) Power spectrum of the force fluctuations. 



there are strong negative events: sometimes the plate is pushed in the opposite direction 
of the mean wind velocity. The spectrum, plotted in figEl^c), is very broad. The thin peak 
at 40 Hz corresponds to the disk rotation frequency. The cut-off at 70 Hz is imposed by 
the filter of the digitizer. Note on the spectrum that the change of slope near 35 Hz is at a 
frequency of the order of U/L where U is the mean velocity at the plate location and L is 
the plate size. 



III. DATA ANALYSIS 



The data analysis of the fluctuations of the normalized force y = F/ Fq has been performed 
in the following way. We first compute the integrated quantity: 

Y{t,T)=- r\{t')dt', (6) 

T Jt 

then we calculate the histograms of the Y{t, r) values in order to obtain the probability 
density functions Tfriy) for each value of the coarse graining time r. The logarithm of the 
functions nriX) are plotted in figO as a function of F for r = Stq, IOtq and 20ro. One 
clearly sees that the PDFs change with the integration time r. For all values of r, they 
are non-Gaussian eventhough the force results from a spatial integration of the pressure 
fiuctuations. The central limit theorem does not necessarily apply here Q| because the 
pressure fiuctuations can be correlated over the plate size: they are caused by turbulent 
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structures of sizes ranging from the integral to the dissipative Kolmogorov scale. In addition, 
one observes the occurence of negative fluctuations, as shown in figI21 
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FIG. 3: Probability density function for the fluctuations of the force on an obstacle placed in the 
von Karman flow, corresponding to the time signals shown in fig|2I^a,b) 

From the PDF 7r^(r), we construct the function Gr{Y) = ln(^^^) for < F < 0.5. 
This small interval around zero allows us to study the symmetry of Gr{Y) even for the 
largest integration times, e.g. t = 20tq. The function Gr{Y) is plotted as a function of Y in 
figEfa) for r = 3tq, 5tq and IOtq. We find that Gr{Y) is very much linear in Y, as predicted 
by eq.l with a slope ^(t) which increases with r. 

We thus turn to the study of the slope a;(r) as a function of r. If eq.l is correct a;(r) 
should be a linear function of r. The slopes ^(t) vs. t are plotted in figEJ^b). We see 
that for t/tq > 4, the dependence of a(r) is linear in r. It is a quite noteworthy feature 
that t/tq = 4 is very close to the rotation period of the disks (this also corresponds to 
the correlation time of the von Karman flow l6j). Indeed, figEJb) shows that the linear 
behaviour of Q:(r), predicted by eq.l, is obtained only for r larger than the mean integral 
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FIG. 4: (a) The function Gr{Y), as defined in the text, plotted as a function of Y; integration 
times: (A) 3tq, (o) 5tq, (□) IOtq. If eqElis verified then Gr(Y) has to be a linear function of Y. 
(b) The function a{T) computed using the values iTt{Y) for —0.5 <Y< 0.5 (o) and 1.5 < y < 2.5 
(*), plotted as a function of r. 

time of the system. This observation will be discussed in more details in the next sections. 

Here we want to stress another aspect of the measurement, namely the importance of 
considering the points around Y = for the study of the Fluctuation Theorem. As a test, 
we have computed the evolution «(t) obtained by the same analysis of the symmetries of 
the PDF around the point Y = 2 instead of the point Y = 0. Corresponding results are 
plotted with * symbols in figlH^b). We observe that in such a case ^(t) is not linear in r. 
We thus find that the point y = has some special property; this is in agreement with the 
FT prediction. 

Summarizing we find that for t/tq > 4 a^r) is linear, and we can write ^(t) = At + B. 
Thus 

ln^^ = (Ar + i?)r (7) 
M-Y) 

We stress that this expression is not in contradiction with eq.l which has been derived for 
r ^ To; rather eqEI should be considered as a non- asymptotic version of eq.l. 

Finally, we check the existence, for the variable of our experiment, of a function indepen- 
dent of T similar to the one defined in eq.3. In order to do this test, we write the PDF as 
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FIG. 5: Sinai scaling function obtained from the PDF measured in the counter rotating disk 
experiment 



T^riY) = At- e"''^^-' Here we use the previously determined function ^(t) instead of r, 
because, as discussed for eq.2, we are not in the asymptotic regime for all values of r. If the 
prediction of ref. Ill "^^^ ^e applied in this case, then C(^) = 1ii[-P(^)Mt-]/«(t) has to be 
independent of r. The functions obtained from the PDFs i^riX) of figOlis plotted in figlHl 
We clearly see that the rescaling is very good. In particular we note that the positive part 
of tIt-{Y), which is strongly non-Gaussian, also scales perfectly. Thus the ^(t) determined 
using eq.3 allows us to rescale the PDF of Y for any r. It is important to stress that the 
rescaling is not possible if r, instead of a(r), is used. Furthermore it is obvious that no scal- 
ing can be observed using the a;(r) obtained by performing the analysis around the point 
Y = 2. 

In order to check the Sinai scaling function we have performed a different data analysis 



as suggested in ref. 



17[. This analysis, which we will not describe here, yields results in 
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close agreement with those plotted in figO 

IV. A WIND TUNNEL EXPERIMENT 

In the previous section we have shown that the PDF of the fluctuations of force on a 
small surface satisfies the symmetries imposed by FT and that they can be rescaled using 
eqOl These two properties were also observed on the PDFs of the local heat flux fluctuations 
in a turbulent thermal convection experiment . It is important to recall that eqEl has 
been proved for global observables and the fact that the symmetry imposed by eqO can be 
observed on a local measurement is not obvious. In addition, the counter rotating disks 
experiment and the thermal convection experiment are spatially bounded flows. Indeed, in 
the von Karman experiment, one observes that, eventhough the fluid is not confined within 
walls, the turbulence remains bounded in a volume whose lateral extent is roughly equal 
to twice the diameter of the driving disks. In such situations, one can argue that a local 
measurement is indicative of the global behavior of the system. 

The question thus arises naturally as whether the same observations can be made in truly 
open systems, where velocity fluctuations are advected outside the observation zone. For 
this reason, we have repeated the experiment on the wind pressure fluctuations inserting 
the measuring surface in the turbulent region of a jet developing inside a wind tunnel. The 
jet is formed by a 5.2 cm nozzle. The obstacle is located 35 diameters downstream of the 
nozzle where the turbulence is fully developed, at distance of 4 diameters off the axis of the 
jet. At the measurement location, the mean wind speed is 1 m/s and the Reynolds number 
based on the jet size is Re = 45000. The force of the wind on the obstacle is measured as 
the bending of the shaft that holds the plate, as detected from the deviation of a laser beam. 
The sampling frequency is I/tq = 100 Hz, and 3.6 10^ data points are accumulated, with a 
16 bit resolution. 

The force fluctuation spectrum is plotted in figEfb). The spectrum is very broad and 
presents no resonance (the peak at 30Hz is due to the wind tunnel driving fan). As in sec. 3 
we define the normalized variable y = F{t)/FQ, and the coarse-grained variable Y after 
integration over a time interval r. The PDFs of Y are plotted in figEI The curves are 
extremely non-Gaussian, with the development of an exponential tail at large values of the 
force, but also important negative excursion on the left-hand side. 
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In order to analyze them, we repeat on the PDFs of the wind pressure in the wind tunnel 
the analysis developed in section 3. Specifically we verify that the previously defined function 
Gr{Y) is linear in Y for —0.5 < Y < 0.5. In figHIa) we plot the measured slope a('r) as 
a function of r. We find that the affine behavior previously observed, ^(t) = At + B, is 
again verified here even for small values of r. It means that the PDFs of the wind tunnel 
measurements do verify the functional form of the fluctuation theorem (eq|21), in an interval 
close to F = 0. 

However we find that for small r values, it is not possible to derive a unique function C(^) 
independent of r which would yield a complete rescaling of the PDFs. This is evidenced in 
figlHIb) where it can be observed that some rescaling property begins to form at the larges 
values of r, e.g. for r > 20tq. It is very interesting to note that this time corresponds to the 
time of flight of the wind from the jet outlet to the measuring point; for times longer than 
this time of flight, the fluctuations can be associated to global variations of the turbulent field 
impacting the obstacle, as a block. This point deserves more investigations, some currently 
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FIG. 6: (a) Schematic diagram of the wind tunnel experiment, (b) and (c) Typical time series of 
the force fluctuations at Re = 45000. (d) Spectrum of the force fluctuations measured in the wind 
tunnel at Re = 45000. 
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FIG. 7: Wind tunnel experiment. PDF of Y for t/tq = 1, 10, 20, 30, 40, 50, upper to lower curve. 




FIG. 8: Wind tunnel experiment, (a) a{T) and (b) C(^)) estimated as in sec. 3 
under progress. 



V. CONCLUDING REMARKS 

In this paper we have described two experiments in which we have measured the pressure 
fluctuations of the wind on a small surface. One of the two experiments, the counter rotating 
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disks, is a spatially correlated flow |16| and the other, the wind tunnel, is an open flow. 
The main results of this investigation are that in the spatially correlated flow the PDFs 
of the fluctuations of force on an obstacle, which are strongly non Gaussian, verify the 
symmetry imposed by the Gallavotti-Cohen fluctuation theorem, and in addition can be 
rescaled according to the Sinai theory. More precisely in order to rescale the PDF one has 
to use the measured «(t) and not simply r as indicated in eqOl The linear function a(r) 
may be seen as a flnite time correction of eqOl which is valid only for r oo. The existence 
of a scaling function C(^) independent of r is very important because it strengthen the 
check of the FT. Indeed a strict test of the FT using eqEl can only be done for values of Y 
close to zero, and one always observe that the probability of getting negative fluctuations 
in an actual system is actually very small. However our measurements show that using the 
value of a(r) estimated for Y close to zero, it is possible to scale all the PDFs onto a single 
function C(^) fo^' This is a non trivial prop erty which was already observed in the 

heat fluctuations in turbulent thermal convection [l^. This property is not observed, at 
least for short times, in the jet data described in sec. 4 of this paper. The main difference 
between the jet experiment, the thermal convection and counter rotating disk experiments is 
that the the flrst behaves as a spatially uncorrelated flow because perturbations are always 
advected downstream of the nozzle. As we have already mentioned this difference is an 
important one. Indeed eqI2]has been proved only for global observables, therefore one can 
picture that if the turbulence is strong enough, a quasi local measurement can provide a 
signiflcant measurement of all of the flow, if the later is spatially correlated. This is not 
necessary true in a open flow where the recurrence time can be very long. This claim is of 
course preliminary and it merits to be studied in more details. 

Another important point which deserves a study is the prefactor in eq|2] which in the 
standard derivation is the phase space contraction rate. In our experiments we have to 
study how this prefactor changes as a function of the Reynolds number. This is needed in 
order to have a more quantitative comparison between theory and experiment. 

Finally we have to stress that the study of the probability of extreme negative fluctuations 
can be very useful in many applications. If one could prove the universality of eqUJin various 
geometries, then one would readily use it to compute the probability of the rare, extremely 
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negative, events. 
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